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First MINOS Results with the NuMI Beam 

B. J. RebeP for the MINOS Collaboration 
a Fermilab, PO Box 500, Batavia, IL 60510, USA 

The results from the initial exposure of the MINOS detectors to neutrinos produced by the Fermilab NuMI 
beam are reported here. The exposure consisted of 1.27 x 10 20 120 GeV protons incident on the NuMI target. The 
data show the observation of 215 neutrinos with energies below 30 GeV while 336 ± 14.4 events were expected. 



The data are consistent with disappearance via neutrino oscillations with IA77123I = 2.74lg'2g x 10 
and sin 2 2<9 23 > 0.87 (68% C.L.). 
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There is substantial evidence that oscillate 
into v T |l|2|3j . The Main Injector Neutrino Os- 
cillation Search (MINOS) was designed to study 
this hypothesis. The probability that a of en- 
ergy E remains a after traveling a distance L 
is 

Pu^u T =l-sm 2 26 23 sm 2 (l.27\Am 2 23 \L/E), (1) 

where jAm^l has units of eV 2 , L has units of 
km and E is in GeV @|. MINOS tests the oscilla- 
tion hypothesis by making two measurements of a 
beam of produced in the Neutrinos at the Main 
Injector (NuMI) beam at Fermilab. The hrst 
measurement occurs at the Near Detector (ND) 
located onsite at Fermilab, 1 km from the pro- 
duction point; the second measurement is made 
at the Far Detector (FD) located 735 km away in 
the Soudan Underground Mine in Soudan, Min- 
nesota, USA. MINOS extracts the oscillation pa- 
rameters by comparing the reconstructed energy 
spectra of the at the ND and FD. 

The NuMI beam is produced using 120 GeV 
protons from the Main Injector. The protons are 
delivered in 10 /is spills, each of which contains 
up to 3.0 x 10 13 protons. Two parabolic horns are 
pulsed with 200 kA of current in order to focus the 
produced ir + and K + toward the detectors. The 
horns are pulsed in such a way as to maximize 
the production of neutrinos in the energy range 
of 1-3 GeV. A total of 1.27 x 10 20 protons on tar- 
get (POT) were delivered for the data described 
below. The neutrino beam is made of 92.9% u^, 
5.8% V a , 1.2% v e and 0.1% V e . 



The MINOS detectors are steel-scintillator 
tracking calorimeters with toroidal magnetic 
fields averaging 1.3 T 5 . The steel planes are 
2.54 cm thick and the scintillator is mounted to 
the steel. The scintillator planes are made of 
4.1 cm wide and 1 cm thick strips. The strips 
in each plane are rotated 45° from the vertical 
and the strips in successive planes are rotated 90° 
from each other. The light produced in the scin- 
tillator is collected in 1.2 mm wavelength shifting 
fibers embedded in the scintillator. The fibers 
transport the light to the multi-anode photomul- 
tiplier tubes (PMTs). The detectors were made 
as similar as possible in order to cancel the major- 
ity of the uncertainties in the neutrino interaction 
modeling and detector response. Both detectors 
yield 6-7 photoelectrons per plane for normally 
incident minimum ionizing particles. The main 
design differences between the two detectors are 
due to the much higher rate (~ 10 5 times) in the 
ND than in the FD. 

The FD is 705 m below the surface, has a 
mass of 5.40 kton, and is composed of 484 in- 
strumented planes. The detector has a regular 
octagonal cross section and is 8 m wide. The scin- 
tillator is read out at both ends of the strips using 
Hamamatsu M16 PMTs. The front end readout 
electronics are designed to provide high precision 
timing information. The FD data were blinded 
until the procedures for event selection and en- 
ergy spectrum prediction were defined and under- 
stood. The blinding procedure hid a substantial 
and unknown fraction of the events in the FD. 
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The ND is 103 m below the surface, has a mass 
of 0.98 kton, and is composed of 282 planes. The 
planes have an irregular octagonal cross section 
and are 4 m tall and 6 m wide. The geometry 
of the planes optimizes containment of hadronic 
showers and allows for the magnetic field to be 
similar to that in the FD. The front end electron- 
ics in the ND are designed to handle the high rate 
of interactions observed with each Main Injector 
spill. Up to 10 interactions can be observed in 
the ND for each spill of protons. As such the 
first step in the reconstruction of the ND data is 
to use timing and spacial information to separate 
the individual interactions. The data in the ND 
were not blinded. 

The energy of each neutrino interaction is 
found in the same way for both detectors. Muon 
tracks are found and their curvature in the mag- 
netic field is fit to determine their energy. The 
hadronic showers are also found and their energy 
is determined. The events selected in both de- 
tectors were required to have visible energy, E V i S , 
less than 30 GeV and the events had to have a 
negatively charged track, a requirement chosen 
to select only interactions. A fiducial volume 
was defined to contain the hadronic energy of the 
event and reject background cosmic ray muons. 
The events were also required to occur within a 
50 /is window surrounding the spill time. The 
background due to cosmic ray muons in the ND 
is negligible. However, an additional constraint to 
reduce the cosmic ray background was imposed in 
the FD; the direction of the reconstructed track 
at its vertex had to be within 53° of the neutrino 
beam direction. The background due to cosmic 
ray events in the FD is estimated to be < 0.5 
events (68% C.L.) as there are no events occur- 
ring within the 50 /is window. 

As the MINOS detectors record both neutral 
current (NC) and CC interactions, a particle 
identification parameter (PID) was defined to de- 
termine whether an event was NC or CC. The 
PID incorporated the probability density func- 
tions for the event length, fraction of energy 
contained in the track and average track pulse 
height per plane to separate the events into these 
two categories. Figure Q] shows the PID for the 
ND and FD data overlaid with the Monte Carlo 
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Figure 1. Data and Monte Carlo distributions 
for the PID parameter in the ND (top) and FD 
(bottom). The dotted vertical lines show the po- 
sition of the cut for CC-like events and the arrows 
indicate the selected events. 



NC and CC distributions. The events with PID 
> —0.2 were categorized as CC events in the FD; 
in the ND the requirement was PID > —0.1. The 
values of the PID were chosen so that the purity 
of the samples in the ND and FD were both about 
98%. The efficiencies for selecting CC events 
with E V i S < 30 GeV in the fiducial volume are 
74% for the FD and 67% for the ND. 

The measured energy spectrum in the ND is 
used to predict the unoscillated spectrum in the 
FD. The method used by MINOS to predict the 
FD spectrum is the Beam Matrix method [5] . In 
this method the ND data are used to measure 
effects such as beam modeling, neutrino interac- 
tions and detector response that are common to 
both detectors. The beam simulation is used to 
derive a transfer matrix that relates in the two 
detectors via their parent hadrons. The matrix 
element gives the relative probability that 
the distribution of secondary hadrons producing 
the observed of energy Ei in the ND will pro- 
duce the observed v„ of energy Ej in the FD. 
The reconstructed energy spectrum in the ND 
is translated into a flux by correcting the sim- 
ulated ND acceptance and then dividing by the 
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Figure 2. Comparison of the FD spectrum with 
predictions for no oscillations for both analysis 
methods. Also shown are the best fit hypoth- 
esis from the Beam Matrix method and the es- 
timated NC background. The last bin contains 
events with 18 < E vis < 30 GeV. 



Figure 3. Confidence intervals for the best fit 
using the Beam Matrix method. The intervals 
include the effects of the systematic uncertainties. 
Also shown are the contours from the previous 
highest precision experiments |l|2!3j . 



calculated cross-sections for each energy bin. The 
flux is multiplied by the matrix to yield the pre- 
dicted unoscillated FD flux. The final step in the 
process is to do the inverse correction for cross- 
section and FD acceptance, resulting in the pre- 
dicted visible energy spectrum. As a cross check 
of this method, the prediction was also done with 
the ND Fit method, which minimizes differences 
between the ND data and Monte Carlo by modi- 
fying the parameters associated with neutrino in- 
teractions and detector response. The FD Monte 
Carlo is adjusted by using the best-fit values of 
those parameters. 

The FD data set contains a total of 215 events 
with E V i s < 30 GeV compared to the unoscillated 
expectation of 336.0 ± 14.4. The uncertainty is 
due to systematic uncertainties associated with 
(a) the fiducial mass calculation and POT count- 
ing accuracy (4%), (b) the hadronic energy scale 
(11%) and (c) the NC component (50%). The ob- 
served energy spectrum is shown in Fig. [2] along 
with the predicted spectra from both methods. 

The data were fit to the hypothesis of — > v T 



oscillations, and the fit has sin 2 26*23 > 0.87 at the 
68% C.L. The best fit for the mass squared differ- 
ence is | Ami, | = 2.741^26 x 10~ 3 eV 2 /c 4 with the 
fit probability of 8.9%. The 68% and 90% confi- 
dence intervals are shown in Fig. [3] along with the 
90% contours from the previous highest precision 
experiments. 
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